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Abstract 


A methodology is presented which uses commercial design and analysis software and the Internet to 
perform concurrent multidisciplinary optimization. The methodology provides a means to develop 
multidisciplinary designs without requiring that all software be accessible from the same local network. 
The procedures ate amenable to design and development teams whose members, expertise and respective 
software are not geographically located together. This methodology facilitates multidisciplinary teams 
working concurrently on a design problem of common interest. Partition of design software to different 
machines allows each constituent software to be used on the machine that provides the most economy and 
efficiency. The methodology is demonstrated on the concurrent design of a spacecraft structure and 
attitude control system. Results are compared to those derived from performing the design with an 
autonomous FORTRAN program. 
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Nomenclature 


B state-space control influence matrix 
Cperf performance output influence matrix 

d t ith structural design variable outer diameter 
d m . x maximum allowable outer diameter 
d mm minimum allowable outer diameter 
D r modal open-loop damping matrix 
e pointing error vector 
E( ) expected value of ( ) 

G f , position gain matrix 
G,. rate gain matrix 
I rxr 1' X r identity matrix 

J objective function 

J r modal inertia matrix 

K r modal stiffness matrix 

L ;j Cholesky factor of position gain matrix 

L, Cholesky factor of rate gain matrix 

M ia3x maximum allowable mass 

M srr structural mass 

P imx average control power 

q r modal amphtude 

t time 

Tr( ) trace of ( ) 
u control input vector 
y p position measurement vector 

y a performance measurement vector 
y r rate measurement vector 
r modal control influence matrix 
O r retained open-loop eigenvector matrix 



1. Introduction 


Multidisciplinary optimization provides a means to integrate two or more design disciplines under the 
same optimization umbrella. The design discourse includes all design parameters concurrently. This 
allows a means of determining the tradeoffs between design parameters in the constituent design 
disciplines. Furthermore, the design parameters are varied until they satisfy a particular performance 
criteria [1-8], Hence, system level performance objectives and constraints can be achieved. This paper 
presents a procedure for performing multidisciplinary optimization using the Internet and UNIX 
commands. The paper illustrates the design procedure on a spacecraft design problem that involves the 
disciplines of attitude control system design and structural design. 

Performing multidisciplinary optimization via the Internet has many advantages. Ref. 7 (Messac et al) 
has demonstrated that multidisciplinary design and optimization has matured beyond the research domain 
to industrial use. Messac has demonstrated a more general approach to multidisciplinary design by 
implementing die design with commercial software packages. The goal of this research effort has been to 
extend multidisciplinary design and optimization such that it is amenable to design and development 
teams whose members and respective expertise and software are not geographically located together. The 
Internet provides a means of communication, transfer of data and transfer of process control. The 
advantage of such an approach is that it allows for more efficiency and economy. For example, this 
approach does not require that one design software package reside on a particular mainframe computer 
simply because the other software is there or putting the more computationally intensive software on a less 
computationally effective machine because other necessary design software is on the machine. Because 
software license costs are dependent upon the machine where the software is installed, software can be 
placed on the platform that will result in the least expensive license. Hence, logistically, this allows for a 
better use of resources. 

The method presented in this paper uses the UNIX operating system [9]. Although not presented, this 
method can similarly produce results with the commercial software packages on computer systems with 
VMS, Windows NT and UNIX operating systems. On all operating systems, it is necessary to use the 
TCP/IP protocol [10]. 

The philosophy used in the development of this design methodology was to extend the existing design and 
analysis methodology beyond implementation using the same local file server (i.e., computers networked 
together such that they share the same files) while making such methods suitable for incorporation into 
modern industry engineering procedures. Although only the facets of control system design and structural 
design are presented here, the procedure used in this paper can be easily be applied to other 
multidisciplinary design and analysis (i.e. thermal, acoustic, plume impingement, etc.) to make a 
combined computer aided design (CAD) and computer aided engineering (CAE) environment. 
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2. Design Optimization Procedure 


The integrated design problem is that of incorporating the two or more disciplines concurrently instead of 
sequentially developing each facet of the system design. All design parameters from each discipline are 
placed in a system level design vector, d . An objective function is prescribed such that as it is reduced in 
value the overall system performance improves. Furthermore, the objective must be explicitly or 
implicitly dependent upon a set of design parameters. Design parameters are given such that they can be 
varied to change the overall performance (design objective). Thus, the design objective is to 


subject to 


Minimize F(d ) 

Gj (a) < 0 j = l,m 


and 

a, < <a u i = l,n 

1 u i 


( 1 ) 

( 2 ) 


where F(CL) is an objective function which, when minimized, will result in a desired performance of the 
system. The design vector, CC, contains n system parameters which are varied through the iterative 
optimization process. Gj (OC) is the jth constraint on the design parameters. There are ifl constraints. 

Each design parameter, <2, , is bounded by lower and upper limits, 0C h and £X U , respectively. System 

physical limitations in hardware or software such as controller torque output (magnitude vs. frequency), 
measurement sampling rate, measurement sampling range (i.e., limits), structural loads limits, structural 
resonance, hardware size, etc. can either be included in the constraint equations or in the design 
parameter upper or lower bounds. 


A flowchart of the optimization procedure is shown in Fig. 1. The procedure is iterative. The first 
iteration requires the user to develop an optimization objective, constraint equations, and bounds for the 
elements of the design vector, hi the first iteration, the objective function is evaluated. If the objective 
value is less than some desired value, the process is complete. If the objective value is greater than some 
desired value, the optimization process searches for new values for the elements of the design vector, OC , 
which will reduce the value of the objective function and satisfy the constraint equations. The permissible 
values for elements of OC are those that are within the limits of the upper and lower bounds of OC . After 
the constraints are satisfied, the design vector is used to generate a new system design. Next, the new 
system produces another response. The response is used to generate another objective function evaluation 
that is then compared to the desired objective value. Iteration of the optimization process is continued 
until some desirable value of Fid) is achieved or the number of iterations has reached some prescribed 
maximum. Finite differencing is used to develop design gradients. 
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3. Design Process Management and Control 


To implement the design process on an internet level, it is necessary to have a means to regulate all 
machines as necessary. The design is regulated by using process status interrogation, copy to remote 
machine and executing remote commands. User privilege (i.e., having a valid account) must exist on all 
machines. Each constituent machine should have a “.rhost” file established which will allow remote shell 
commands to be mutually executed among all machines. This is done by placing the Internet protocol 
(IP) address and user name for each system used. An example of information which needs to be placed in 
the “.rhost” file is given below: 


(IP address A) (username on A) 
(IP address B) (username on B) 


(IP address N) (username on N) 


The example given above considers the case where each constituent piece of software needed for each 
discipline (e.g., A, B, ... N) resides on a different machine. The entire process is regulated by a single 
FORTRAN or C program, hr the example which will be presented, all programs are written in 
FORTRAN and all operating systems are UNIX. 

hr the first iteration: 

1. The regulating program reads a tile containing the initial design vector. 

2. If the design objective is not achieved, the optimization program is called to generate new 
design parameters. 

3. The new design vector is partitioned to files which are used as inputs to the design software at 
each respective internet site. 

4. Design and analysis software is executed at each site. 

5. Objectives and constraints are evaluated using results from the design software. 


In succeeding iterations, step 2 through 6 are repeated until the optimization is complete or the 
optimization has gone through a predetermine number of iterations (Fig. 1). The regulating program has 
a system level call to the UNIX procedural files. These procedure files execute design and analysis 
software. It is through the procedure files that the analysis required in the optimization can be 
multiplexed to different machines at remote sites. 

For example, the complete system is compose of two disciplines, A and B (e.g., control design and 
structural design). All software for discipline A resides on the machine A. All engineering and technical 
expertise for discipline A is located at site A. Similarly, all software and expertise for discipline B resides 
at site B. The optimization software resides at A. The logistics for implementing analysis on machine A 
is relatively simple. All design parameters used in A are saved to a file, the analysis software for A is 
then executed. Results from the analysis are saved into a file. Similarly, the design parameters used in B 
are saved to a file, design_var_B. When the optimization regulating program requires analysis results 
from site B, a system call inside the program such as: 
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1 

2 42 

3 

4 

5 

6 

7 

8 

9 

10 
12 

13 

14 

15 31 


call system! ‘get-resultS-from-B’) 
call system! “ read_status_B ’ ) 
open (3 ,file='job_status' ) 
read(3,*) jobstat 
close(3) 

ifljobstat.eq. 1.0) tlien 
call system! ‘read_results_B’) 
go to 3 1 
else 

call system('sleep 10') 
go to 42 
endif 
continue 


The procedural file ‘get_results_from_B' contains the following instructions 
echo “0.0” > 'jobstatus’ 

rep “design var B” “122. 133.45.67:design_var_B” 
rsh 122. 133.45.67 run B software 


The instructions above execute a remote process, run_B_software, on the machine B with EP address 
122. 133.45.67. The procedural file run_B_software exists on B and does the following: 

echo “0.0” > 'job_status' 

designsoftwareB 

wait 

echo “1.0” > 'job_status' 

The design and analysis software at site B, design_software_B, uses the design var iables that were placed 
in the file, design_var_B, runs its analysis and places results in the file results_B. The “wait” command 
prevents other processes from running until the command preceding the “wait” is finished. The 
optimization regulating program Hines 2- 14) will constantly interrogate the file job_status using the 
procedural file read_status_B which contains the following command: 

rep "122.133.45.67:job_status" job status 

If the contents of job_status is 0.0, the program sleeps (i.e., remains idle) for 10 sec. and then re- 
examines the file job_status. If the contents of job_status is 1.0 the optimization regulating program 
copies the results from machine B to machine A using the procedural file read_results_B which contains 
the following instructions: 

rep "122.133.45.67:results_B” "result_B" 

The results from the analysis process from machines B and A can be used to evaluate constraints and 
objectives and the optimization process continues. 
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4. Numerical Example: Integrated Design of the Earth Pointing System 


4.1 System Modelling 

The methodology of partitioning the design problem to various Internet remote sites has been 
demonstrated using the system presented in Ref. 1 (Maghami et al). The results used in Ref. 1 are also 
compared to those produced in this paper. Results produced by Maghami et al used an aggregate 
FORTRAN program for structural design, control design and optimization. The system used in Ref. 1. 
and this paper is the NASA Earth Pointing System (EPS). The EPS is shown in Fig. 2. The EPS is a 
generic geostationary platform that was used as a means to compare the results of various control and 
structural design methodologies. Although the model is fictitious, it is a realistic representation of 
existing Earth observing spacecraft. The model also has coupling between the structure and the control 
system. Any change in inertia distribution will effect the controller response inversely. If the controller 
bandwidth is higher than the fundamental (lowest) structural frequency, then there could be possible 
adverse interaction between the controller and the structure. A concurrent multidisciplinary optimization 
approach will allow designers to take advantage of interaction. The EPS system is designed for fine 
attitude control and vibration suppression. It is a 30-m-long truss structure consisting of 10 structural 
bays, a 7.5 m diameter radial antenna and a 15 m diameter radial antenna. All structural members, 
including the antenna, are hollow tubes with circular cross section and 1.59 mm thickness. 

The Automated Design Synthesis (ADS) (Ref. 11) program is used for optimization. All control analysis 
is evaluated using Matlab (Ref. 12). NASTRAN (Ref. 13) is used for all structural analysis A FORTRAN 
program regulates the optimization process. Two machines with different IP addresses are used in this 
demonstration. NASTRAN, ADS and the regulating program reside at one of the addresses. The control 
software on Matlab resides at the other address. 


Finite element modeling using NASTRAN is used to generate the structural model in physical space. 
Usually the models are very large. Modal truncation is used to reduced the size of the model. Modes 
beyond the controller bandwidth are eliminated in the modal model. Ref. 1 presents a more detailed 
discussion of model formulation and development of controller design. The system equations in modal 
space using modal coordinates are 


J r q r + D r q r + K r q r = Oj Bu = T u 


(l) 


where q r is a r X 1 vector of modal amplitudes; J r , D r , K r , are, respectively, the generalized 
inertia, damping, and stiffness matrices; O r is a r X n matrix whose columns are the r open-loop 

eigenvectors associated with the included modes; B is a n X m control influence matrix and U is the m 
X 1 control input vector. If mode shapes are normalized with respect to the inertia matrix and modal 
damping is assumed then J r = I rxr , D r = Diag 2( ), K r = Diag 

(&r,&r,...,fy~), where (O t and ^ ( / = l,2,...,r) are the open-loop frequencies and damping ratios, 
respectively. 


=Tq r 

y r = r q r 

[q 

y pert ^ ' peif 


( 2 ) 
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The position and rate measurements used in the feedback loop are_y ; , and y r , respectively. The part of 
the output used in the performance evaluation is given by_y pelf . 

u = "G p y p - G,y, 

G P = KK 
g, = l,i; 

(3) 

The static dissipative (Ref. 1) controller position and rate gain matrices are G p and G ; . , respectively. 
The elements of their respective Cholesky factorization matrices, L p and L ; . are used as control design 
variables. 


4.2 EPS Optimization 

Minimization of the steady-state RMS pointing error at the large antenna is used as a design objective. 
Because geostationary satellites must maintain very high pointing tolerances, the objective is a relevant 
design concern. The objective is represented by the following function: 

/ = lim^rr^?[c(f)c / (0]j) 

(4) 

Average control power , P irax , limitation must be considered in real-world applications and can be 
included as a design constraint 


jiin(Tr{£[u(()u r (()]})< ^ 

(5) 

Tr is a trace of its argument and E is the expected value of its argument. The average control 
power, P mK , is limited to 3N 2 - m 2 . Launch vehicle payload is also a design concern. Spacecraft mass 
must be less than the maximum allowable for the launch vehicle. Furthermore, any savings in spacecraft 
mass can be used to increase the amount of fuel on the spacecraft and thus give the spacecraft a longer 
life. Therefore structural mass is included as constraint 


M s , < M 


max 


( 6 ) 

Although structural truss diameters can be used as design variables, their size must be within some 
bounds to maintain the structural integrity of the spacecraft and to also bound the spacecraft size. 
Therefore 


(L.. < d. < d. 


(7) 
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5. Discussion of Results 


The method presented in this paper is compared to the design produced in Ref. 1 by Maghami et al. Ref. 
1 uses a single FORTRAN program for concurrent optimization of the EPS control system and structure. 
This paper produced results using separate software packages for control design, structural design and 
optimization. The regulation of the optimization process was done via Internet commands. The Internet 
also was used a means of passing data between processes. The final spacecraft total mass in both designs 
had a 1.4% difference, hi this study, the mass was reduced from 5760.1 Kg to 5676.0 Kg. Both designs 
final average control power differed by only 0.07 N 2 -m 2 . Ref. 1 produced a RMS pointing error that was 
2.1 rad. larger than that produced in this paper. This method reduced the pointing error from 71.55 
microradians to 14.68 microradians. The attitude gain Cholesky matrix produced using both methods 
differed by as high as 16.2%. However, the elements for the attitude rate gain Cholesky matrix were 
within 5.6% of each other. The major differences between the two methods were in the final structural 
design. 

Both designs resulted in stronger (larger diameter) large antenna mounts and weaker small antenna 
mounts. The design produced in this paper had smaller truss diameters in section 3 than Ref. 1. The 
smaller diameter resulted in more structural flexibility near the small antenna, hi section 2, only the 
longhorn diameters differed significantly. The section 1 longhorns produced by Magma et al were twice 
the diameter produced in this paper. The section 1 battens produced by the Internet design 1.8 times those 
of Ref. 1. 

hi Figs. 3 through 11, the optimization design history is presented. The EPS RMS pointing error (Fig. 3) 
converges to the final value of 14.68 rad. after approximately 240 iteration steps. During the course of the 
design the mass and average control power constraints (Figs. 4 and 5) are violated between steps 75 and 
270. The rate and attitude gam elements follow the same trends. At approximately step 100 and 210, 
they increased their values and the values remain relatively constant. The structural variables follow the 
same trend (Figs. 8-10) until iteration step 150. After step 260 they converge to their final values. The 
small antenna mount diameters (Fig. 1 1) reach their final value at approximately step 97. After step 260, 
process converged on a final large antenna mount design (Fig. 1 1). The process is continued beyond step 
270 to allow all constraints to be satisfied. 

The total time involved in the parceling of the total design to different machines and the design 
orchestration via the Internet was very insignificant compared to the run-time of the structural design 
software (NASTRAN) 
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6. Concluding Remarks 


A method has been presented that demonstrates how multidisciplinary optimization can be accomplished 
using two or more computers with design discipline partition among the computers using the Internet. 
Use of the Internet provides a means for multidisciplinary design teams to interact remotely to accomplish 
a system level design. The system level design process was accomplished using a FORTRAN program 
that regulated the design analyses on all the computers involved. The method presented in this paper used 
the UNIX operating system. The TCP/IP protocols were used for data transfer and remote computer 
process control. Each discipline design was accomplished using commercial software. 

Structural and control design of the Earth Pointing System (EPS) was used to demonstrate the 
methodology. The system is fictitious but possesses the detail of a realistic Earth observing spacecraft. 
Results from the design process were compared to those produced by Magma et al (Ref. 1). Ref. 1 used 
single FORTRAN program for control design, structural design and multidisciplinary optimization. The 
final design produced using the method presented in this paper. Average power control, total mass and 
structural member size were included as design constraints. The design objective was to minimize the 
RMS pointing error of one the antenna. 

The final design resulted in a better pointing performance than Ref. 1 . The final power required and the 
final total spacecraft was within 1% of that produced by Ref. 1. Other final design results were 
commensurate with Maghami et al. This method has demonstrated that multidisciplinary design can be 
automated and regulated among different computers. The results produced also show the viability of using 
concurrent multidisciplinary design partitioned to remote internet sites. 
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Appendix Summary of Procedure Files, Matlab M-files and Fortran Programs 
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Fig. 2 Earth Pointing System. 
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Fig. 6 Diagonal Elements of Cholesky Matrix of the attitude gain matrix. 
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Fig. 11 Large and small antenna mount diameters. 
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